Quarkonium formation time in quark-gluon plasma 
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The quarkonium formation time in a quark-gluon plasma (QGP) is determined from the space- 
time correlator of heavy quark vector currents using the quarkonium in-medium mass and wave 
function obtained from heavy quark potentials extracted from the lattice QCD. It is found that 
the formation time of a quarkonium increases with the temperature of the QGP and diverges near 
its dissociation temperature. Also, the quarkonium formation time is longer if the heavy quark 
potential is taken to be the free energy from lattice calculations for a heavy quark pair, compared 
to that based on the more negative internal energy. 
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I. INTRODUCTION 



Quarkonia of heavy quark-antiquark bound states are 
promising probes for the QGP created in relativistic 
heavy- ion collisions Many studies have been carried 
out in recent years. These include experimental measure- 
ments of its yield in heavy ion collisions at both the Rel- 
ativistic Heavy Ion Collider (RHIC) (H-Q ^^id the Large 
Hadron Collider (LHC) as well as theoretical stud- 

ies based on various models These studies have 

indicated that a quantitative description of quarkonium 
production in these collisions requires the understand- 
ing of its interactions in both the produced QGP and 
the initial cold nuclear matter. One of the important 
QGP effects is the dissociation of quarkonium by thermal 
patrons [13, IHj . The resulting thermal dissociation 
width of a quarkonium is determined by its dissociation 
cross section and the density of partons. In both the 
leading order (22| - [2^ and the next-to-leading order cal- 
culations [25| in perturbative Quantum Chromodynam- 
ics (pQCD), the quarkonium dissociation cross section 
is characterized by its dipole nature through the deriva- 
tive of the quarkonium wave function with respect to the 
relative momentum between the heavy quark and anti- 
quark. However, the thermal dissociation of quarkonium 
in QGP is relevant only after the quarkonium is com- 
pletely formed from the heavy quark-antiquark pair that 
are produced from initial hard scatterings between collid- 
ing nucleons. While the heavy quark pair are produced 
during the time of 1 /tuq with toq being the heavy quark 
mass, a quarkonium takes a longer time to be formed 
from the pair. This is because the typical time scale 
needed to separate the bound state from the continuum 
or the excited states are related to the inverse of their en- 
ergy differences, which in the heavy quark limit of QCD 
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scales as the bound state energy mgg^, where g is the 
QCD coupling constant [2^ . 

There have been several atteinpts to determine the for- 
mation time of quarkonium |26l428l | . One of them is based 
on the space-time correlator of the heavy quark current 
operator This approach has the advantage that all 
parameters involved in the calculation can be related 
to experimentally measured quantities in heavy quark 
production from electron-positron annihilation, e"''e^ — > 
QQ. In the present study, we extend this method to 
study the quarkonium formation time in QGP by us- 
ing the quarkonium in-medium mass and wave function 
that are determined by taking the heavy quark poten- 
tial to be either the free energy or the internal energy 
of a heavy quark pair in QGP as extracted from lattice 
calculations j29l - l3ll |. 

This paper is organized as follows: In SecHIl we briefly 
review the quarkonium formation time in the approach 
based on the space-time correlator of the heavy quark 
vector current. We then describe in Sec. Illll the quarko- 
nium formation time in QGP by using the quarkonium 
in-medium mass and wave function calculated from the 
lattice heavy quark potentials. Finally a summary is 
given in Sec. IIVI 



II. FORMATION TIME OF QUARKONUM 

For the heavy quark vector current operator J^{x) 
QjiiQ, its space-time correlator H^,y(x) is given by 



(27r)4 

Using the dispersion relation 
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to relate the real part of the heavy quark correlator in 
energy-momentum space or pair polarization function 
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n(g2) to its imaginary part, one can rewrite, for x 7^ 0, 

Eq. (p as m,!!!,!!! 
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where 
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is the relativistic causal propagator in the coordinate 
space with Ki being the modified Bessel function and 
T being the EucHdean proper time [33 |. 

The imaginary part of the heavy quark pair polariza- 
tion function is related to the ratio of the cross section 
for heavy quark pair production to that for dimuon pro- 
duction in electron-positron annihilation, 



Imn(s) = 
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QQ,s) 
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The cross section ratio R(s) = '^M 1 in the 

above equation can be parameterized by the sum of res- 
onances and a continuum 12811: 
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where Ji, Mi, T1 ^ , and Ti are the spin, mass, dielectron 
and total decay widths of resonance i, respectively; eg is 
the electric charge of the heavy quark and the threshold 
energy ^sth is taken to be twice the mass of open heavy 
flavor. 

Substituting Eqs. ^ and (O into Eq. ([3]), one ob- 
tains following contributions to the heavy quark space- 
time correlator from resonances and the continuum 1281: 
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These expressions show that the contribution from the 
continuum dominates at earlier times and is then gradu- 
ally taken over by that from resonances. This is consis- 
tent with the picture that the QQ pair produced in e^e~ 
annihilation is initially formed in a continuum state and 
then evolves into a quarkonium state. 

By taking the ratio of Fi{T) = ni(T)/n>(T), where 
ni(T) is the heavy quark space-time correlator due to 
resonance i and n> (r) is the total heavy quark pair po- 
larization function excluding contributions from quarko- 
nium resonances with masses lower than the mass of 



quarkonian i, one can then obtain the distribution of the 
formation time of quarkonium i according to [28j 
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III. 



QUARKONIUM FORMATION TIME IN 
QGP 



Although the quantities Mi, , and ^sth appear- 

ing in Eq. ^ are well known if one is interested in deter- 
mining the formation time of a quarkonium in free space, 
their values in hot dense medium depends on theoret- 
ical models. In the present study, we determine them 
using the heavy quark potential taken from the poten- 
tial energy between a heavy quark-antiquark pair that 
is extracted from lattice calculations |29l - l3ll |. While the 
heavy quark potential based on their free energy is di- 
rectly extracted from lattice calculations, that based on 
their internal energy depends also on the entropy density. 
As a result, the internal energy potential is stronger than 
the free energy potential. At present, it is still an open 
question as to which potential energy should be used in 
determining the properties of a quarkonium in QGP. We 
therefore use both potential energies for the present study 
and compare the resulting predictions. 

For a given potential between a heavy quark and anti- 
quark pair, the Schrodinger equation has the form 



2mQ-— V^-f y(r,r) 



^,{r,T) = M,^,{r,T), (9) 



where mg is the bare mass of heavy quark Q and 
^i (r, T) is the wave function of quarkonium i at tem- 
perature T . Introducing the potential energy at infin- 
ity, V{r = 00, T), the Schrodinger equation is modified 

to Hi! 



1 



—V'^ + V{r,T) 
rriQ 



= -l^2mQ + V{r ^ 00, T) - M,^^,{r,T) 



= -^2mQ-M,}^,ir,T) 
= -eiil;.i{r,T), 



(10) 



where V{r,T) = V{r,T) - V{r = 00, T) and it vanishes 
at infinity, ffiQ = mq + V{r ^ oo,T')/2, and £i is the 
binding energy of quarkonium i at temperature T . 

From its wave function at the origin, the dielectron 
decay width of quarkonium i can be calculated according 
to 111 



M? 



|^(r = 0)p. 
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where a is the fine structure constant. 
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FIG. 1: (Color online) Wave functions of T (IS) and T (2S) 
for the case of free energy potential at T = 1.0, 2.0, 2.4 Tc 
and at r = 1.0, 1.05, 1.08 Tc, respectively. 



FIG. 2: (Color online) Distribution of T (IS) formation time 
for the case of free energy potential at T = 1.0, 1.4, 1.8, and 
2.2 Tc. 



As an example, we consider the formation time of bot- 
tomonia in QGP. Their in-medium masses and wave func- 
tions are obtained by solving Eq. ([TU)) with the bottom 
quark mass of rub = 4.65 GeV. Figure [T] shows their wave 
functions of T(IS') and T{2S) for the free energy poten- 
tial at several temperatures. It is seen that as tempera- 
ture increases, the wave function of quarkonium becomes 
broader. As a result, the probability for the heavy quark 
and antiquark in the quarkonium to annihilate and decay 
into dielectron becomes smaller. We note that the bind- 
ing energies of T(IS') and T{2S) are larger and their wave 
functions also have larger values at the origin if the more 
attractive heavy quark potential based on the internal 
energy is used. 

Taking the threshold energy for the continuum contri- 
bution in Eq.(I7]) to be -y/s^ = Smg, we have calculated 
the distribution of the T (IS) formation time at several 
temperatures for the case of free energy potential, and 
they are shown in Fig. [51 The distribution peaks at small 
r at T = 1.0 Tc but spreads to later times at high tem- 
perature. There are two reasons for this behavior of the 
T formation time in QGP. First, the smaller F^^*^ at 
high temperature delays the time for the contribution of 
quarkonium to the heavy quark correlator n(T) to be- 
come dominant. Secondly, the small binding energy of 
quarkonium at high temperature makes its formation to 
take longer time. In the large Euclidean time r limit, the 
contributions to the heavy quark correlator from a res- 
onance i and the continuum are proportional to e"*^'"^ 
and e"^'^'^, respectively [l^. Because the resonance 
mass Mi is always smaller than the threshold for the 
continuum, the contribution from a resonance eventually 
becomes dominant. However, if the binding energy of 
quarkonium is very small, it takes a long time for the 
contribution of the quarkonium state to become much 



more important than that of the continuum states. 
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FIG. 3: (Color online) Average formation times of T(15), 
T(2S'), and T(3S') as functions of temperature for the cases 
of free energy (upper) and internal energy (lower) potentials. 

Figure [3] shows the average formation times of T(15), 
T(25'), and T(35') as functions of temperature, 



(Tfo 



SdT tP,{t) 
JdrP^ir) 



(12) 



for the two cases of free energy and internal energy poten- 
tials. It is seen that the average formation time of quarko- 
nium increases with temperature and diverges near its 
dissociation temperature. This behavior is similar to that 
of the quarkonium radius in QGP [l^. This is reason- 
able because as the quarkonium radius increases at high 
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temperature, it takes longer for the quarkonium to be 
formed (26j . Since the quarkonium is less bound and has 
a larger radius in the case of the free energy potential 
than that of the internal energy potential, the formation 
time is thus longer in the former than in the latter case. 
Compared to the formation time of 0.32 fm/c obtained in 
Ref. for a T(IS') in free space, our values of 0.44 fm/c 
and 0.3 fm/c at T = Tc for the cases of free energy and 
internal energy potentials are slightly larger and similar, 
respectively. 



1.0 




T/T 
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FIG. 4: (Color online) Ratio of the mean distance between 
bottom quark and antibottom quark in T (IS) to its formation 
time as a function of temperature for the cases of free energy 
(solid line) and internal energy (dashed line) potentials. Solid 
squares are the temperature dependence of the QCD coupling 
constant (scaled by 1.1) at the screening distance between 
heavy quark and antiquark pair extracted from the lattice 
free energy. 

In Fig. m we show the ratio of the mean distance be- 
tween bottom quark and antibottom quark in T (IS) to 
its formation time as a function of temperature, calcu- 
lated according to / d'^xr|i/'(?')p/(Tform)- It can be in- 
terpreted as the average relative velocity between bot- 
tom and antibottom quarks before they form the T (IS) 
state. Results in Fig. |4] show that the relative velocity is 
smaller for the case of free energy potential than for that 
of internal energy potential, and that it decreases as tem- 
perature increases. The latter is consistent with the fact 
that the wave function of bottomonium in momentum 
space has a small relative momentum at high tempera- 
ture [36j . However, this does not necessarily mean that 
only bottom quark pairs with small relative momentum 
contribute to bottomonium formation at high tempera- 
ture, because even if a heavy quark pair initially have 
a large relative momentum, they can lose energy dur- 
ing propagation through QGP [33, [13 and thus become 
slowly moving relative to each other. 

It is useful to discuss the result shown in Fig. |4] using 
the relevant scales for the heavy quark system in QCD. 



The heavy quark system is characterized by three scales; 
hard toq, soft p = rngv and ultra-soft E = rnqv^ in 
terms of the heavy quark mass and velocity. The different 
scales are the basis for the potential non-relativistic QCD 
(pNRQCD) approach, which is an effective field theory 
for heavy quarks with bound states (39| . In the mg — > 00 
limit, one notes that the soft and ultra-soft scales are 
proportional to the inverse Bohr radius 1/ao and the 
binding energy of the system, which scale with the QCD 
coupling as rnqg^ and rnqg^, respectively [22l |. There- 
fore, the heavy quark pair mean distance to the quarko- 
nium formation time ratio shown in F ig. [H should be pro- 
portional to or aes{T). In Ref. [io], aqq{rscrecn,T) 
was extracted from the lattice free energy by estimating 
the nonperturbative coupling constant at the screening 
distance between the quark and antiquark pair, which 
should be the most relevant distance for the heavy bound 
state. It is quite interesting to see that the temperature 
dependence of Q;gg(rscrocn, T), scaled by an overall factor 
of 1.1, is close to the ratio in Fig. 2] extracted from the 
free energy potential. 



IV. SUMMARY 

We have studied the quarkonium formation time in 
QGP by using the approach based on the space-time 
correlator of heavy quark vector currents. The imagi- 
nary part of the resulting heavy quark pair polarization 
function, which is the spectral function of heavy quark 
pair in e"^e~ annihilation, is constructed by solving the 
Schrodinger equation with the heavy quark potential ex- 
tracted from lattice calculations. The real part of the 
polarization function, which is related to its imaginary 
part by the dispersion relation, then provides the infor- 
mation on how different states of the heavy quark pair 
evolve with time. 

Using bottomonia as examples, we have found that the 
average formation time of a quarkonium from a heavy 
quark-antiquark pair increases with temperature and di- 
verges near the dissociation temperature of the quarko- 
nium. Furthermore, the quarkonium formation time is 
longer if the heavy quark potential is taken to be their 
free energy instead of their internal energy from the lat- 
tice calculations. We have also found that the average 
relative velocity between the heavy quark pair before 
they form the quarkonium, calculated via the ratio of the 
mean distance between the heavy quark and antiquark 
in a quarkonium to its formation time, decreases with 
increasing temperature. Our results thus indicate that 
because of the longer formation time of quarkonium at 
high temperature, to fully describe quarkonium produc- 
tion in relativistic heavy ion collisions requires a proper 
treatment of the medium effects on heavy quark pairs in 
the QGP before they form the quarkonium bound state. 
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